Model-independent source imaging using two-pion correlations 
in 2 to 8^4 GeV Au + Au collisions 
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We report a particle source imaging analysis based on two- 
pion correlations in high multiplicity Au + Au collisions at 
beam energies between 2 and 8^4 GeV. We apply the imaging 
technique introduced by Brown and Danielewicz, which al- 
lows a model-independent extraction of source functions with 
useful accuracy out to relative pion separations of about 20 
fm. The extracted source functions have Gaussian shapes. 
Values of source functions at zero separation are almost con- 
stant across the energy range under study. Imaging results 
are found to be consistent with conventional source parame- 
ters obtained from a multidimensional HBT analysis. 



One of the robust predictions of Quantum Chromody- 
namics is that strongly interacting matter can exist in 
a state with colored degrees of freedom, i.e., quarks and 
gluons, if subjected to sufficiently high temperature or 
density ffl. This deconfined state of matter, often called 
Quark-Gluon Plasma (QGP), may have been present in 
the Universe a few microseconds after the Big Bang. It 
is also believed that QGP can be produced in collisions 
of heavy ions at high energies. Experimental searches 
for evidence of QGP are prominent and challenging pri- 
orities in modern heavy-ion physics. A major difficulty 
is that the sought-after QGP is a transient state which 
persists only for timescales on the order of fm/c. After 



expansion and cooling, any QGP necessarily undergoes 
a phase transition to hadronic matter, possibly leaving 
little evidence behind in the final state of the collision. 
Numerous experimental observables have been proposed 
as signatures of QGP creation in heavy- ion collisions [||. 
One of the more promising of these is based on the ex- 
pectation that the larger number of degrees of freedom 
associated with the deconfined state should manifest it- 
self in increased system entropy, which ought to survive 
subsequent hadronization and freeze-out. The specific 
prediction is that QGP formation will cause the system 
to expand more and/or interact for a longer time and pro- 
duce more particles, relative to a system which remained 
in the hadronic phase throughout the collision process. 

Intensity interferomctry pi has been extensively used 
to extract information about spatial-temporal properties 
of heavy-ion collisions. The HBT approach ||]-[|] is well 
understood from the theoretical point of view, and is the 
mainstream analysis method for meson (ir, K) intensity 
interferometry. Meson final-state interactions are sim- 
ple, so one may remove most of the distortion of the 
final-state wavefunctions with a simple Coulomb correc- 
tion. Furthermore, in the traditional analysis one does 
not ask for more than the widths of the underlying source 
function. A novel approach — the imaging technique 
introduced by Brown and Danielewicz [|| O] - offers 
the opportunity to use any class of like-pair correlations 
to reconstruct the entire source function, as opposed to 
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the conventional model-dependent parametrization. The 
two approaches have natural connections, which are es- 
pecially clear for the case of non-interacting spin-zero 
bosons In this paper, we report the first application 
of the source imaging technique to experimental data, 
specifically negative pions from Au + Au collisions at 
beam energies between 2 and 8A GeV, with particular 
emphasis on verification and study of the expected con- 
nection with the source parameters obtained from HBT 
analysis. Comparing the two methods not only allows us 
to test the model underlying the HBT parametrization, 
but also lends credibility to the source imaging approach 
in other contexts such as proton-proton correlations. 

We present data from experiment E895 H], in which 
Au beams of kinetic energy 1.85, 3.9, 5.9 and 7.9A GeV 
from the AGS accelerator at Brookhavcn were incident on 
a fixed Au target. The analyzed 7r~ samples come from 
the main E895 subsystem — the EOS Time Projection 
Chamber (TPC) ||l5|| , located in a dipole magnetic field 
of 0.75 or 1.0 Tesla. Results of the E895 ir~ HBT analysis 
were reported earlier fl6| , and in this Letter, we focus on 
an application of imaging to the same datasets. 

In order to obtain the two-pion correlation function 
Ci experimentally, the standard event-mixing technique 
was used 17 1. Negative pions were detected and recon- 
structed over a substantial fraction of 47r solid angle in 
the center-of-mass frame, and simultaneous measurement 
of particle momentum and specific ionization in the TPC 
gas helped to separate ir~ from other negatively charge 
particles, such as electrons, K - and antiprotons. Con- 
tamination from these species is estimated to be under 
5%, and is even less at the lower beam energies and higher 
transverse momenta. Momentum resolution in the region 
of correlation measurements is better than 3%. Event 
centrality selection was based on the multiplicity of re- 
constructed charged particles. In the present analysis, 
events were selected with a multiplicity corresponding to 
the upper 11% of the inelastic cross section for Au + 
Au collisions. Only pions with pt between 100 MeV/c 
and 300 MeV/c and within ±0.35 units from midrapid- 
ity were used. Another requirement was for each used 
7r~ track to point back to the primary event vertex with 
a distance of closest approach (DCA) less than 2.5 cm. 
This cut removes most pions originating from weak de- 
cays of long-lived particles, e.g. A and K°. Monte Carlo 



simulations based on the RQMD model 18 indicate that 
decay daughters are present at a level that varies from 
5% at 4 A GeV to 10% at 8 A GeV, and they lie preferen- 
tially at pt < 100 MeV/c. Finally, in order to overcome 
effects of track merging, a cut on spatial separation of 
two tracks was imposed. For pairs from both "true" and 
"background" distributions, the separation between two 
tracks was required to be greater than 4.5 cm over a 
distance of 18 cm in the beam direction. This cut also 
suppresses effects of track splitting [ pi] . 

Besides such cuts, the correlation functions themselves 



were corrected for Coulomb final-state interactions and fi- 
nite momentum resolution effects, using an iterative tech- 
nique similar to that used by the NA44 Collaboration 
[19]. See [161 for details. 

Imaging |Q|JlCf1 was used to extract the shape of the 
source from the measured correlation functions. It has 
been shown p| that imaging allows robust reconstruc- 
tion of the source function for systems with non-zero life- 
time, even in the presence of strong space-momentum 
correlations. The main features of the method are out- 
lined below; see Refs P, ^0| , p^[ for more details. The full 
two-particle correlation function may be expressed as: 



C P (Q)-1 = J drK(Q,r)S P (r). 



(1) 



where K — \^^\r)\ 2 — l, the total momentum of the pair 
is P, the relative momentum is Q — pi — P2, and the 
relative separation of emission points is r. Both r and Q 
are written in the pair center of mass frame. Here, <frq 
is the relative wavefunction of the pair and for spin-0 
bosons in the absence of final-state interactions, it is the 
symetrized sum of free waves: 
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However, in general it may be quite complicated. The 
goal of imaging is to the determine the relative source 
function (Sp(r) in Eq. (0)), given Cp(Q). The prob- 
lem of imaging then becomes the problem of inverting 
K(q, r). The source function, Sp(r), is the distribution 
of relative separations of emission points for the two par- 
ticles in their center-of-mass frame. 

In this letter, we analyze the pion correlations in Qmv 
The transition from the full three-dimensional problem 
in Eq. (|l|) to the angle-averaged problem is straightfor- 
ward: one expands the source and correlation in spherical 
harmonics (|] and keep £ ~ m — terms. The angle- 
averaged version of Eq. (^) is 



C(Q lnv ) -1 = 4* drr 2 K(Q mv ,r)S(r) 



(3) 



Here Qi nv = \J Q 2 — Qo and since Qq = in the pair 
CM frame, Q ntt , = |Q|. Thus, while the source one re- 
constructs is a distribution in the pair CM frame, the in- 
version itself may be done with correlations constructed 
in any frame. In this equation, the kernel is simply aver- 
aged over the angle between q and r: 



K(Qi 
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d(cos6 QT )K(Q,r) 



(4) 



For identical spin-zero bosons with no FSI, this kernel is 



{.Qinv ? " Sin \Qinv'^)/Qinv^' t 



(5) 



2 



We comment that since the r = point of the source 
is a maximum, it is spherically symmetric in a neigh- 
borhood around that point and S(r — ► 0) = S(r — > 0). 
We now proceed as in plj and expand the radial de- 
pendence of the source in a Basis Spline basis: 5(r) = 
X)j S 3 Bj{r). With this, Eq. (g) becomes a matrix equa- 
tion d 



KijSj with a new kernel: 



(6) 



Imaging reduces to finding the set of source coeffi- 
cients, Sj, that minimize the x 2 ■ Here, x 2 = Yli (@i ~ 
KijSj) 2 1 A 2 Ci. This set of source coefficients is 



- T,il( KT ( A2c )~ lK )~ lKTB ]ji( C i - 1) where K T 
is the transpose of the kernel matrix. The uncertainty 
of the source is the square-root of the diagonal ele- 
ments of the covariance matrix of the source, A 2 S = 
(K T (A 2 C)~ 1 K)~ 1 . 
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FIG. 1. Measured two-pion correlation functions for Au + 
Au central collisions at four beam energies. Lines represent 
Gaussian fits to the data. 



For the case of noninteracting spin-zero bosons with 
Gaussian correlations, there is a natural connection be- 
tween the imaged sources and "standard" HBT parame- 
ters jl2| — the source is a Gaussian with the "standard" 
HBT parameters as radii: 



[R 2 } = 
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Eq. (^) is the most general Gaussian one may use, but 
usually one assumes a cylindrical symmetry of the single 
particle source so that there is only one non-vanishing 
non-diagonal element R 2 ol . The asymptotic value of the 
relative source function at zero separation S(r — > 0) is re- 
lated to the inverse effective volume of particle emission, 
and has units of fm -3 : 



S(r - 0) 



A 



(2V5F)Vdet [i? 2 ] 



(9) 



As it was shown in [jl2| S(r — * 0) is an important 
parameter needed to extract the space-averaged phase- 
space density. Figure ^ shows measured angle-averaged 
two-pion correlation functions for central Au + Au col- 
lisions at 2, 4, 6 and 8A GeV. Figure || shows relative 
source functions S(r) obtained by applying the imag- 
ing technique to the measured two-pion correlation func- 
tions. Note that the plotted points are for representation 
of the continuous source function and hence are not sta- 
tistically independent of each other as the source func- 
tions are expanded in Basis Splines |2l| . Since the source 
covariance matrix is not diagonal, the coefficients of the 
Basis spline expansion are also not independent which is 
taken into account during x 2 calculations. The sources 
shown in Figure || are consistent with zero, within errors, 
in the region past ~ 20 fm. Thus with this technique, we 
may reconstruct the distribution of relative pion separa- 
tions with useful accuracy out to r ~ 20 fm. 
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Here A is a fit parameter traditionally called the chaotic- 
ity or coherence factor in HBT analyses and [R 2 ] is the 
real symmetric matrix of radius parameters 



r, fm 

FIG. 2. Relative source functions extracted from the pion 
correlation data at 2, 4, 6 and 8A GeV. 
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The images obtained at each of the four E895 beam 
energies are rather similar in shape, and upon fitting 
with a Gaussian function, values of x 2 per degree of free- 
dom between 0.9 and 1.2 are obtained. Results of fits to 
the source function and correlation functions are shown 
in Table |. One can see that source radii extracted via 
both techniques are similar, further confirming the valid- 
ity of a Gaussian source hypothesis. Indeed, the some- 
what smaller errors on the Rg parameters compared to 
Rc2 are consistent with the Gaussian fit function being 
in closer overall agreement with the data points in the 
case of the imaging approach. 

TABLE I. Radius parameters of Gaussian fits to the 
extracted source (Rs) functions and measured correlation 
functions(i?c2) for different energies. 



E b (AGeV) 


2 


4 


6 


8 


fts(fm) 


6.70±0.04 


6.35±0.03 


5.56±0.03 


5.53±0.05 


R C 2 (fm) 


6.39±0.2 


6.05±0.1 


5.51±0.15 


5.61±0.28 



Figure ^| compares the effective volumes of pion emis- 
sion inferred from the standard Bertsch-Pratt pion HBT 
parametrization (open circles) with the effective volumes 
derived from the image source functions shown in Fig. ^ 
(solid circles), recalling that S(r — ► 0) = S(r — > 0). Es- 
sentially, this figure plots the inverse of the left- and right- 
hand sides of Eq. ||. It can be seen from Fig. || that the 
agreement between imaging and the HBT parametriza- 
tion is fairly good. Values of source functions at zero 
separation as well as HBT source size parameters (pub- 
lished in Ref. jl6| and reproduced in Table [n]) are approx- 
imately constant within errors across the 2 to 8A GeV 
beam energy range. 

TABLE II. Fit parameters for the Bertsch-Pratt HBT pa- 
rameterization of the pion correlation functions for E895 beam 
energies used for evaluation of S(r — > 0). 



E b (AGeV) 



6 



A 

Ro (fm) 
R s (fm) 
Ri (fm) 
Rli (fm) 



0.99T0.06 0.74±0.03 0.65T0.03 0.65T0.05 

6.22±0.26 5.79±0.16 5.76±0.23 5.49±0.31 

6.28±0.20 5.37±0.11 5.05±0.12 4.83±0.21 

5.15T0.19 5.15±0.14 4.72±0.18 4.64±0.24 

-2.43±1.71 0.43T1.03 2.17±1.20 -0.65il.85 



In summary, we present measurements of one- 
dimensional correlation functions for negative pions emit- 
ted at mid-rapidity from central Au + Au collisions 
at 2, 4, 6 and 8A GeV. These correlation functions 
are analyzed using the imaging technique of Brown and 
Daniclcwicz. It is found that relative source functions 
S(r) have rather similar shapes, while zero-separation 
intercepts S(r — ► 0), which are related to the effective 
volume of pion emission, are almost constant across the 
range of bombarding energies under study. Distributions 
of relative separation have been measured out to 20 fm, 
and the extracted source functions are approximately 
Gaussian. We have performed the first experimental 
check of the predicted connection between imaging and 
traditional meson interferometry techniques and found 
that the two methods arc in good agreement. Overall, 
results of this investigation strengthen the rationale for 
applying systematic imaging analyses to a variety of par- 
ticle species (i.e., strongly interacting particles such as 
protons, antiprotons, etc., as well as pions), and open the 
door to more detailed inferences about spacetime struc- 
ture. Such inferences can play an important role in re- 
solving outstanding questions about deconfined nuclear 
matter. 
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FIG. 3. Values of inverse relative source functions at zero 
separation as a function of beam energy, obtained using imag- 
ing (solid circles) and HBT (open circles). 
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